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This paper discusses the modes of optical resonators, and optical modes 
of propagation or Gaussian beams of light. The passage of Gaussian beams 
through lenses, telescopes, sequences of lenses, and lenslike media is studied. 
Mode matching formulae are derived. A complex beam parameter is intro- 
duced for which the law of transformation by any given optical structure 
can be written in the simple form of a bilinear transformation (A BCD law). 
Resonators with internal optical elements and their transmission line duals 
are also investigated. Effective Fresnel numbers and curvature parameters 
are determined, which allow one to infer the diffraction losses, the resonant 
conditions, and the mode patterns of the various systems. Results are ob- 
tained for resonators with internal lenses, sequences of lenses with irises 
inserted between the lenses, resonators with internal lenslike media, trans- 
mission lines consisting of a lenslike medium with periodically spaced 
irises, and resonators with one very large mirror. 

I. INTRODUCTION 

The theory of Fresnel diffraction is the basis for an understanding of 
optical resonators 1 " 5 and of optical modes of propagation. 23,4 Fresnel 
diffraction explains the mode patterns and diffraction losses of optical 
resonators, and the beam waist and spreading of the modes of propaga- 
tion or "Gaussian beams." In this paper we will discuss how these 
Gaussian beams of light are transformed on their passage through lenses, 
telescopes, various lens combinations, and lenslike (guiding) media, 
and how these optical systems affect the properties of optical resonators 
when inserted between the resonator mirrors. 

We will assume that no additional aperture diffraction effects are 
introduced by these optical systems, i.e., that the apertures of the 
internal lenses can be regarded as infinitely large. The imaging laws of 
geometrical optics are therefore expected to apply, and we will use them 
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wherever possible, as they generally simplify the algebraie derivations 
and at the same time provide some physical insight. 

Some of the problems to be investigated here in greater detail have 
already been treated in the literature. Goubau 8 has given some mathe- 
matical relations between the parameters of Gaussian beams trans- 
formed by a thin lens. The recently published mode matching formulae 7 
are the result of a computation which will now be presented. Resonators 
with internal lenses have also been discussed in the literature, 8-11 and 
we have used the concept of an effective distance 9 - 10 in a previous publica- 
tion. 9 In several cases an alternative to our algebraic approach is the 
graphical method of Collins, 11 who introduced the circle diagram 1112 
for Gaussian beams. 

In the following we will first establish the rules of imaging for Fresnel 
diffraction with attention to the imaging of the phase fronts which are 
of particular importance for optical modes. Then we will list expressions 
for the focal length and the principal planes of various optical systems of 
interest, because these parameters are needed later for application of the 
imaging rules. This listing includes the parameters of the telescope, of 
sequences of lenses, and of sections of lenslike medium. Armed with 
these tools we will study the passage of Gaussian beams through lenses 
and various optical systems. The paper is concluded by an investigation 
of optical resonators with internal optical elements and their transmis- 
sion line duals. Effective Fresnel numbers and curvature parameters 
arc determined which allow one to infer the diffraction losses, the reso- 
nant conditions, and the mode patterns of the various systems. Results 
are obtained for resonators with internal lenses, sequences of lenses with 
irises inserted between the lenses, resonators with internal lenslike 
media, transmission lines consisting of a guiding medium with periodi- 
cally spaced irises, and resonators with one very large mirror. 

II. IMAGING RULES 

While geometrical optics deals with rays, the theory of Fresnel dif- 
fraction deals with (scalar) fields. To describe the field distribution, we 
use complex amplitudes E(x,y,z) and a Cartesian (x,y,z) coordinate 
system. We consider a wave that propagates in the direction of the 
optic axis (z axis). Within the assumptions of Fresnel diffraction an 
ideal thin lens of focal length / transforms the incoming wave with a 
field E Mt (x,y,z = const) immediately to the left of the lens into a 
wave with the field 

~J k ' 2/ ' / ^ 
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immediately to the right of the lens. Here k is the propagation constant. 
The thin lens produces a phase shift which is proportional to the square 
of the distance to the optic axis, while the intensity distribution is the 
same on both sides of the lens. 

Consideration of spherical waves provides a link between (1) and 
the laws of geometrical optics according to which a spherical wave with 
a radius of curvature 2?i at the left of the lens is transformed into a 
wave with curvature radius 7? 2 as shown in Fig. 1 . The radii i?i and R 2 
are related by 



(1/Zfc) + (1/B.) = 1//. 



(2) 



For Fresnel diffraction the transverse field distribution of a wave with 
a spherical phase front of radius R is given 2,3 by 



E(x,y,z = const) = exp (—jkr 2 /2R) 



where 



= x* + y\ 



(3) 



(4) 



and R is counted positive for a phase front that is concave if observed 
from the left. For spherical phase fronts of radius R\ on the left and — R 2 
on the right of the lens (where the phase front curvature is negative, 
as shown in Fig. 1) we can express E\ e t t and bright with the help of (3), 
compare the exponents in (1), and find the same relation (2) between 
Ri , R 2 , and / as for the spherical waves of geometrical optics. 

To discuss imaging consider an object, i.e., the field E\(x\ , yi) in an 
object plane, and its image E 2 (x 2 ,y 2 ) in the corresponding image plane 
(see Fig. 2). The distances di and d 2 between the lens and the two 
planes are related by 



(lM) + (l/rf.) = l//. 



(5) 



, — PHASE FRONTS 




Fig. 1 — Lens transforming phase front of spherical wave. 
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IMAGE PLANE 

*2 




Fig. 2 — Imaging of field distribution by a thin lens. 

We know from geometrical optics that the intensity distributions of 
the object and the image are similar. This is, of course, still true for 
Fresnel diffraction by any field aperture in the object plane. Assuming 
that no aperture diffraction effects are introduced by the thin lens, one 
can use the Fresnel diffraction formula to relate Ei and E 2 (see Appendix 
A ) and arrive at 



E 2 (x 2 ,y 2 ) = -%E X 



( d x di \ 



exp 



- jk ( 



' /,+ * + |£ 



D 



(6) 



with r 2 ? ~= x 2 + y£. The factor d x /d 2 in this equation follows from 
conservation of energy; the arguments — (d\/d 2 )x 2 and — (di/d 2 )ij 2 
indicate that the image is inverted and magnified by d 2 /di . The first 
.^two terms in the exponent are simply due to the phase shift k(di + d 2 ) 
which the light wave suffers in propagating from the object to the image 
plane, while the third and last term is of particular importance for our 
considerations. It describes an additional phase shift proportional to 
r 2 which appears in the field distribution of the image. Apart from 
this additional phase shift the amplitude and phase distribution of the 
image and the object are scales of each other. 

The expression for the additional phase shift follows also from geo- 
metrical optics (see e.g. Appendix B), and it is related to the thick- 
mirror formulae, 13 as we shall see later. It is also obtained by studying 
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the passage of Gaussian beams through a lens. The additional phase 
shift does not appear in Abbe's theory of imaging; he finds that the 
image is strictly similar to the object, both as regards the amplitude 
and phase distribution. 14 But Abbe used the Fraunhofer diffraction 
theory, where phase terms proportional to r are neglected. 

For Fresnel diffraction the r dependence of the additional phase 
shift suggests that one should use spherical reference surfaces instead 
of plane ones, as shown in Fig. 3. By proper choice of the curvature of 
these surfaces tangential to the image and object planes, one can achieve 
an image field on one surface that strictly reproduces the object on the 
other surface in amplitude and phase. For an object reference surface of 
radius R\ and an image reference surface of radius Ri one gets for the 
fields additional phase factors of exp (~jkri 2 /2Ri) and exp {—jkr?/2Ri), 
respectively. These phase factors cancel the additional phase shift in 
(6) if 



1 

'/7 a 



}_<tf , ldi 

Rid* 2 f do/ 



(7) 



After some algebraic manipulations involving (5) this relation can be 
rewritten as 



1 



+ 



1 



1 



(8) 



dt + B, (h - Ri 
This simply means that the center of curvature C\ of the object surface 
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Fig. 3 — Imaging of fields with spherical wave fronts; centers of curvature 
are images of each other. The corresponding spherical reference surfaces are used 
when fields with nonsDherical phase fronts are imaged. 
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is imaged onto the center of curvature C 2 of the image surface. Thus, 
whenever the centers of curvature of the image and object surfaces are 
images of each other we have an image which is a strict (scaled) reproduc- 
tion of the object as regards both the amplitude and phase distribution, 
with no additional phase shift. 

The imaging rules discussed above can also be used to study imaging 
by a combination of lenses (or by any optical system that can be re- 
garded as such). It is not necessary to apply the rules step by step to 
each individual thin lens of the combination. It is generally simpler to 
determine the parameters of the equivalent thick lens as usual in geo- 
metrical optics. The place of / is then taken by the combined focal 
length of the system, and object and image distances (di and d 2 ) are 
measured from the principal planes of the thick lens. 

III. FOCAL LENGTHS OF VARIOUS OPTICAL SYSTEMS 



3.1 The Ray Matrix 

When one traces a paraxial ray through combinations of lenses and 
lenslike media, the quantities of interest are the position Xi and the 
slope xi of the ray in the input plane, and the corresponding quantities 
x 2 and x 2 in the output plane (see Fig. 4). There is in general a linear 
relation 181617 between the output and input quantities which can be 
written in matrix form as 



(9) 



We will call this A BCD matrix the ray matrix of the system. Because 
of reciprocity the determinant of the ray matrix is generally unity: 

AD - BC = 1. (10) 

It is easy to determine the focal length and the principal planes from 
the elements of the ray matrix of an optical system. By tracing a beam 
that leaves the output plane parallel to the optic axis (x 2 ' = 0) we find 
the location of the focal point on the input side. Its distance Si from 
the input plane is obtained as 



x 2 

1 

x 2 


= 


A B 
C D 




Xi 

i 

Xi 



Si = — -, 

Xi 



i 2 '=o 



7) 
C 



(11) 



where we refer to Fig. 4. Similarly, we find for the distance s 2 between 
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Fig. 4 — Reference planes for optical system. 
the output plane and the corresponding focal point 



x 2 



x, '=(1 



4 



(12) 



To find the principal plane on the input side we follow an input ray 
from the focal point until its distance from the axis is equal to the 
position x-i of the corresponding output ray and have 



h = 



•c 2 — Xl 



lo '=o 



(13) 



where the distance h { between the principal plane and the input plane 
is measured positive as shown in Fig. 4. On the output side we find 
similarly 



h* = 



Xj — Xi 

a- 2 ' 



x,'=(l 



(14) 



The focal length / of the system is obtained by calculating the distance 
between a principal plane and the corresponding focal point 



/ = s t + hi = so + h 2 = 



x{ 



El 



(15) 



Using the linear relations of (9 ) together with the last three expressions, 
one finally gets 



/ = - (I/O 



(16) 



462 



THE BELL SYSTEM TECHNICAL JOURNAL, MARCH 1965 



/>, = (D - \)/C (17) 

h = (A - 1)/C (18) 

where the thick-lens parameters are expressed in terms of the elements 
of the ray matrix. For later use we also write down the matrix elements 
as functions of the lens parameters which follow from the last expres- 
sions 



A = 1 - (V/) 

B = hi + h, - (Mt//) 

C = -(1//) 

D = 1 - (h/J). 



(19) 
(20) 
(21) 

(22) 



3.2 The Two-Lens Combination — Telescope 

The lens parameters of a combination of two lenses are well known 
and are listed here for completeness and for later use. The combination is 
shown in Fig. 5. For lenses of focal lengths /i and / 2 spaced at a distance 
d we have 

l/f = (i/M + (1//1) - WfM 

df 



h 
where the lens planes are used as input and output planes. 



(23) 
(24) 



(25) 



h 






v 

L-h,- 



Fig. 5 — The two-lens combination. 
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For a slightly misad justed telescope the lens spacing is 

d=h+ f, - M 



(20) 



where Ad measures the misadjustment. The lens parameters of the 
telescope can be written as 



,Af2 

3 Ad 



hi = 



fjd 
Ad 



Ad 



(27) 
(28) 
(29) 



3.3 Sequence of Lenses 

A periodic sequence of lenses of equal focal length / and lens spacing 
d is shown in Fig. 6. The reference planes are chosen just to the right 
of each lens. The elements of the ray matrix S of one section of the 
sequence (i.e., one lens spaced at a distance d from the input plane) are 
well known 15 ' 17 and are given by 



a- 



1 



1 - 



(30) 
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Fig. 6 — Sequence of lenses of equal focal length. 
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They relate the position and slope (xi and Xi) of the ray just after the 
first lens to the ray position and slope just after the zeroth lens 



(31) 



The ray to the right of the nth lens is related to the input ray by the 
nth power of the ray matrix of one section 



Xl 


= S 


x 


1 




1 


Xi 




Xo 



X n 




.To 


1 
X n 


= S n 


i 

x 



(32) 



The matrix elements of S" can be computed with the help of Sylvester's 
theorem 18 and are well known. 15,16 One has* 



>- 



1 



sin 

sin nd — sin (n — 1 

— T sin nd 
h 



d sin nd 
— - ) sin nd — sin (n — 1 



0-jD- 



(33) 



with 



cos e = 1 - (rf/2/„). 



(34) 



We can now employ (16) and obtain for the focal length / of n sec- 
tions of a periodic sequence of lenses 



/ = / (sin 0/sin nd). 



(35) 



The distance of the two principal planes from the input and output 
planes (zeroth and nth lens) follows also from (33) with the help of 
(17) and (18). One finds 



and 



hi = id/2) +/(1 - cosnfl), 



h 2 = -(d/2) +/(1 - cos nd). 



(36) 



(37) 



If we measure the distance h of the principal planes from the midplanes 
between the lenses as shown in Fig. 6 we have 



h = /(l — cos nd). 



(38) 



* These matrix elements can be written in terms of Chebyshev polynomials 
of the second kind of the variable [1 - (d/2/o)]. 
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A more complicated sequence of lenses is shown in Fig. 7. Here a 
lens of focal length f\ is followed by a lens of focal length / 2 and vice 
versa. The lens spacings are d] and r/ 2 in sequence as shown in the figure. 
This sequence of lenses can be reduced to the simpler type discussed 
above. We can regard it as a sequence of thick lenses formed by lens 
pairs of focal lengths f\ and f-> . The focal length / ( , of the thick lens is, 
according to (23), given by 



V/n = (l//i) + (1//0 - (rfi//,/,), 



(39) 



and expressions for the principal planes are given in (24) and (25). 
The distance d between the output principal plane of a thick lens and 
the input principal plane of the consecutive thick lens is obtained as 



,/ - d, !- /,, r h-i - d, 4-./W/ 1 (L + I 



(40) 



With the principal planes as reference planes, rays passing through this 
sequence of thick lenses behave the same way as rays passing through 
a sequence of lenses of equal focal length that are equally spaced. We 
can therefore use the expressions (33) and (34) obtained above. With 
(39) and (40) the latter becomes 



cos 0=1 — 



di + d 2 



(W) 



+ 



did* 
2fi/. 



(41) 



3.4 Lenslike Medium 

A lenslike medium or "guiding medium" is one whose refractive index 
n varies near the optic axis as in 



<- — d 2 — *4< d, »j* — d 2 — *ir*- 



— d,— 



-d 2 -*| 



h, 



ift ffi T-f. n 



14' 






Fig. 7 — Sequence of lenses of alternating focal length with alternating lens 
spacings. 
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(—a 



n = wo i - i ,-, 



(42) 



where n is a constant, r is the distance from the optic axis, and b meas- 
ures the degree of the variation of n. A medium of this kind can be pro- 
duced by inhomogeneities in laser crystals ' or by a radial variation 
of the gain in high-gain gaseous lasers. 21 Another important example 
is the medium of the recently reported gas lens. 

To trace rays in a lenslike medium one uses the differential equation 
for light rays. 25 For paraxial rays this ray equation has the form 



no 



dz 2 



— -n = -4n * r „ 
dx o 1 



(43) 



for the distance x(z) of the ray from the z axis. A corresponding rela- 
tion holds for y(z). The solution is, again, a linear relation between the 
ray position and slope in the output plane (x and x ) and the cor- 
responding input quantities Xo and x 



cos 2 r 
b 



b . n z 
2 sm2 b 



2 . z n z 

— T sin 2 - cos 2 r 
/> b b 



•I'd 



Xfl 



(44) 



A typical ray path is shown in Fig. 8. To calculate the optical parameters 
for a section of lenslike medium immersed in a medium with a refractive 
index of unity (vacuum), we invoke Snell's law to relate the ray slopes 



n = i 




4— h— | 



INPUT PRINCIPAL PLANES OUTPUT 

Fig. 8 — Ray path in lenslike medium. 
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at the section boundaries. For paraxial rays we have approximately 

#vac = nax ; .T 0v ac = Wo . (45) 

Now we use (16) and find for the focal length of a section of length I 

b 



/- 



~ ~TT (46) 

2no sin 2 r 
o 



(for n = 1 this formula has been given in Refs. 23 and 26, for example). 
The distance h of the principal planes from the input and output planes 
respectively (see Fig. 8) is computed with the help of (17) and (18). 
One obtains 

h = A tan I . (47) 

2n b 

The above expressions have been derived for a focusing medium where 
6^0. For a defocusing medium we have b 2 ^ 0, and the expression 
for the focal length becomes 

| b | 

* o • u o l ' (48) 

2wn smh 2 r-j-f 

I '' I 
The location of the corresponding principal planes is described by 

h = jAJ tanh * . (49) 

IV. OPTICAL TRANSFORMATION OF GAUSSIAN BEAMS 

4.1 Light Propagation in Free Space 

Near the optic axis an optical mode of propagation or Gaussian beam 
is regarded as a TEM wave with a spherical phase-front and a trans- 
verse field distribution that is described by Laguerre-Gaussian 2 or 
Hermite-Gaussian 3 functions. The two beam parameters of interest are 
the "spot size" or beam radius w(z) and the radius of the phase front 
R(z). In any beam cross section of a fundamental mode the field varies 
as 



^(-^-^m) 
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and is specified by w and R. The light beam expands as it propagates 

9 ^ fi 9fi 97 

through space as shown in Fig. 9. The law of expansion is ' ' ' ' 



2 9 

W = Wo 2 



L \7rw7 J ' 



(50) 



Here the z is measured from the beam waist where the phase front is 
plane and the beam reaches its minimum radius wo . For R(z) we 

J 8,8,6,86,37 

have 



Dividing (50) by (51) we find 



(51) 



(52) 



TW _ XZ 
\R TTWo 2 

which we can use to rewrite the terms in the round brackets, and ex- 
press wo and 2 in terms of w and R 



Wo = / 2\2 

R 


\irw-/ 




(53) 
(54) 



U — PHASE FRONT 



Fig. 9 — Contour of Gaussian beam of light. 
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4.2 Beam Transformation by a Lens 

When a Gaussian beam passes through a lens a new beam waist is 
formed, and the parameters in the expansion laws are changed. Assume 
the light beam is propagating to the right. Before passing through the 
lens it has a beam waist a distance d\ away from the lens with a beam 
radius W\ as shown in Fig. 10. The lens produces another beam waist a 
distance d-> away with a beam radius w 2 . The distances d x and d 2 are 
measured positive as shown in the figure (for a negative di one has a 
virtual waist). In the following we will establish some relationships 
between beam parameters of the incoming beam (identified by the 
subscript 1 ) and the parameters of the transformed beam (subscript 2 ) . 

The far field angles 27 0i and 8 2 of the two beams are computed from 
(50) as 

0i = \/itW\ ; 02 = X/7rt02 • (55) 

From these two angles follow immediately the beam radii w if and w 2f 
in the two focal planes of the lens where the image of the far field appears 

w if — /02 = \J/ttW2 (56-) 



w-i/ = J6\ = \f/irwx . 



(57) 



The beam radius in one of the focal planes is, of course, independent of 
the spacing between the lens and the beam waist of the other beam. It 
follows from (51) that the center of curvature of the far field phase 
front is in the beam waist. According to the imaging rules of Section II, 
corresponding centers of curvature are images of each other (where we 
take the phase fronts as reference surfaces). We therefore have to 
determine the image of a beam waist to find the curvature center of 
the phase front in the focal plane on the other side of the lens. The 
distance d 2 between the lens and the image of the waist w 2 follows from 




Fig. 10 — Gaussian beam transformed by a lens. 
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(1/4') + (1/4) - 1//, (58) 

and the radius of curvature R\ f in the left focal plane is equal to the 
spacing between that image and the focal plane 

Rif = d 2 ' - f. (59) 

Combining (58) and (59) we have 

f 

d, - j 

and correspondingly 



Rif = t J — i (60) 



R v = -t£— 4 (61) 

cli — J 

for the radius of curvature in the right focal plane. Rif and Ry are 
independent of the beam radii w 2 and w 1 , respectively, a fact that can 
be used for mode matching into confocal resonators. 

To relate the beam waists we use (56) and (60) to write 

irwi/ \(d 2 — f) 



XRif TW2 2 

and similarly 

irw 2 / _ X(rfi — /) 
\R 2 f irtoi 2 



(62) 



(63) 



To express w 2 in terms of Wi and d\ we insert (57) and (63) into (53) 
and find 



w 2 

r. 



This relation, first given by Goubau, relates the beam radius of the 
waist of the transformed beam to the parameters of the incoming beam. 
A corresponding relation for the spacing d 2 between the lens and the 
beam waist w 2 is found by inserting (61) and (63) in (54). The result is 

do. - f = (di - /) ~ 7—57, . (65) 



(rfi - /)" + 



(¥) 



The above expressions were derived with the help of the imaging 
rules of Section II. As mentioned before, these rules apply not only to 
thin lenses but also to thick lenses and lens combinations. Therefore, 
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if d\ and do are measured from the principal planes the results given 
above are valid for the transformation of Gaussian beams by thick 
lenses. 

4.3 Mode Matching 

In experiments with optical modes one often wants to transform a 
beam with a given beam radius iv\ at the waist into another beam of 
waist radius w 2 . One wants to "match" the modes of one optical system 
(like a laser resonator) to the modes of another one (an optical trans- 
mission line for example). This can be done by selecting a suitable lens 
and by properly adjusting the waist spacings di and d 2 , where we refer 
to Fig. 10. The proper spacings are given by the mode matching formulae 7 
derived below. 

We combine (62) or (63) with (52) and obtain 

2 



di — f _ wi_ 

(k — f I02 2 

This is used to rewrite (64) in the form 



(66) 



i s --L(4-/)(*-/)+ife*Y. (67) 



•uh 



Multiplying (67) by w*j we arrive at 

(di -f)(d,-f) =f-f* (68) 

where we have defined 

/ = w(wiW 2 /\). (69) 

To arrive at the mode-matching formulae we multiply or divide (68) 
by (66), extract the square root, and find 



u>i 



or 



di - / = ± - 1 V/ 2 - /o 2 (70) 

w 2 



-f=±-\fF=U' (71) 

Wi 



As discussed in Ref. 7, one achieves mode matching by choosing a 
lens (or lens combination) with a focal length / that is larger than / 
or equal to it. For a given lens there are generally two ways open to 
match the modes. One can choose either the plus sign in both (70) and 
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(71 ), or the minus sign. For/ = f there is only one set of proper spacings 
<h = di = / = /o . 

4.4 Complex Beam Parameter — A BCD Law 

In the foregoing we have used two parameters to characterize a 
Gaussian beam in a given beam cross section: the spot size or beam 
radius w, and the radius of phase front curvature R. We define now a 
more abstract complex beam parameter q 

\/q= (1/R) -j(\/rw 2 ). (72) 

The propagation and transformation laws for this beam parameter are 
particularly simple and allow one to trace Gaussian beams through 
more complicated optical structures. The old parameters R and w can, 
of course, be recovered from the real and imaginary parts of 1/q. Note 
that we can regard the circle diagram of Collins 11 as plotted in the com- 
plex plane of the variable j/q, and the circle diagram of Li as plotted 
in the complex plane of jq . 

In terms of the complex beam parameter the laws of propagation 
(50) and (51) have the simple and compact formf 

q = qo + z (73) 

as one easily verifies by inserting (50) and (51) into (72). Here 

go = i(W/X) (74) 

is the complex beam parameter at the beam waist. Because of the 
linearity of (73) the parameters q t and q* of two arbitrary beam cross 
sections are related by 

q-i = q\ + d (75) 

where d is the distance between the two planes of interest measured 
positive in the direction of the optic axis. 

The beam parameters qi and q% to the left and to the right of a lens 
are related by 

l/<? 2 = (l/«i) - (1//) (76) 

which simultaneously states the transformation of the phase fronts as 
in (1) or (2), and the fact that the beam radii (widths) are the same 
on both sides of the lens [compare (1)]. 

t Similar propagation laws for optical modes have been used independently by 
]). A. Kleinman, A. Ashkin, and G. D. Boyd in an analysis of second-harmonic 
generation in crystals and by G. A. Deschamps and P. E. Mast in their recent 
paper in Proc. Symp. Quasi-Optics, Polytechnic Inst. Brooklyn, 1964, p. 379. 
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The imaging law (6) applied to Gaussian beams takes the form 

qi d-2 2 q\ f d 2 

if written in terms of the complex parameters gi and q 2 of the beam in 
the object or image planes, respectively. Comparing with (7) and (8) 
one can also write this relation between the parameters of the object 
and the image as 

+ ^—=7. (78) 



di + qi do — q 2 f ' 

Using (75) and (76) one can easily determine how an incoming beam 
with the parameter qi at a distance di from a lens is transformed. The 
parameter q 2 of the transformed beam at a distance d 2 from the lens is 
obtained as 



g 2 = 



(l " f ) q l + (d 1 + d 2 -f) 



-f + (.-f) 

To establish a link to the transformation laws for the real parameters 
developed before, Ave multiply both sides of (79) with the denominator 
of the right side. Then we postulate that we have beam waists at d\ 
and d 2 by putting q\ = jirWi/\ and q 2 = jirw 2 /\. If we compare the real 
parts of the resulting expression, Ave obtain relation (68), and comparing 
the imaginary parts we find (66 ) . 

Let us now regard qi and q 2 as the beam parameters in the input and 
output planes of an optical system described by its ray (ABCD) matrix 
as in Section III. This system is also described by its focal length and 
its principal planes as calculated from (16), (17), and (18). To relate 
qi and q 2 Ave use (79) and put c/i = h x , and d 2 = h 2 . Comparing Avith 
(19), (20), (21), and (22) we see that 

which AA r e shall call the ABCD law. The q parameters of the input and 
the output are related by this bilinear transformation. The ABCD 
law says that the constants of this transformation are equal to the 
elements of the ray matrix. The ray matrices of several optical struc- 
tures are given in Section III, and Ave shall use the ABCD law to study 
the passage of Gaussian beams through some of these structures. 

There appears to be a very close connection between Gaussian light 
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beams and the spherical waves of geometrical optics. In fact, all the 
important laws of this chapter are formally the same for a spherical wave 
with a radius of curvature q. One is therefore tempted to regard a Gaus- 
sian beam as a spherical wave with a complex radius of curvature. For 
the limit of infinitely small wavelengths the curvature radius becomes 
real and one has a spherical wave of geometrical optics. 

The ABCD law allows also a kind of "black box" approach to the 
study of optical modes. One can, for example, inquire about the mode 
parameters of a sequence of equal black boxes, i.e., optical structures 
characterized by their ray matrix elements A, B, C, and D. For a mode 
the beam parameter at the output of a black box is equal to the param- 
eter at the input (qi = q 2 = q). From (80) follows a quadratic equation 
for the mode parameter q 

Cq + (D - A)q- B = 0. (81) 

The solution of this equation can be written as 

r^^sW 4 -^^ (82) 

from which one can obtain the beam radius or spot size of the mode 
and the radius of curvature of its phase front. 

4.5 Beam Transformation by a Telescope 

In this chapter we shall study the passage of a Gaussian beam through 
a telescope consisting of two lenses of focal length /i and f 2 , respectively, 
spaced at a distance d = /i + fa — Ad. The "misadjustment" Ad is 
assumed small. The focal length and the location of the principal planes 
of the telescope are given in (27), (28), and (29). We consider an in- 
coming beam with a beam radius wi at its waist, and the waist spaced 
at a distance si from the first lens as shown in Fig. 11. We want to de- 
termine the location s 2 of the waist of the outgoing beam and its beam 
radius w 2 . 

The distances of the waists to the corresponding principal planes are 

di = Sl + h ■ d 2 = s 2 + h . (83) 

From this we find with (24) and (27) 

* l =$ + ¥(* -l). (84) 

Inserting this expression together with (27) in (64) we get for the beam 
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waist w-i 



Fig. 11 — Gaussian beam passing through a telescope. 
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:s5: 



which is correct to first order in Ad. We see that the ratio of the beam 
waists is more or less equal to the ratio of the focal lengths of the lenses. 
There is only a slight dependence on the position of the input beam waist 
for Ad 7* 0. 

To determine the location of the output beam waist we use (84) 
and the corresponding expression for (d^/f) — 1 to rearrange (08) as 



Sl — fa . 6'2 — fa 



fa 2 



w 



Ad 
/i/t 



[6-06") 



+ 



hk 



¥fi- 



(8c; 



Inserting (85) and expanding to first order in Ad this becomes 

Ad[ , .,2 /W\H 



— /i , s 2 — fa 

i* ft ' 



fa 2 



(87) 



For a well-adjusted telescope we have Ad = 0, and the distances 
between the beam waists and the focal planes of the corresponding 
lenses (i.e., si — fa and s 2 — fa) scale like the squares of the focal lengths 
of the two lenses. The signs in (87) indicate that for an input waist 
which lies to the left of the focal plane of the input lens one has an out- 
put beam waist to the left of the focal plane of the output lens, and 
conversely for an input waist to the right of the input focal plane. 
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4.6 Beam Transformation by a Sequence of Lenses 

Consider now a sequence of lenses of equal focal length / spaced at a 
distance d as shown in Fig. 6. Immediately to the right of each lens an 
optical mode of this structure has a phase front with a radius of curva- 
ture of —2/ and a beam radius w m given 3 ' 4 ' 26 by 



* l JU-l = ^ (88) 



ttwJ 2/ V d d 

where is denned in (34). To the right of each lens the complex beam 
parameter (72 ) of a mode is therefore 

I * t -j±!. (89) 

q m 2/ d 

Assume that a Gaussian beam is injected into the lens sequence, and 
call its complex beam parameter in the input plane </i . If gi = q m , 
then we have launched a mode of the system, and the parameter of the 
beam to the right of every lens is q m . For qi 9± q m we use the A BCD law 
(80) to compute the beam parameter g 2 to the right of the nth lens. The 
elements of the ray matrix of n sections of the lens sequence are given 
in (33), and we use them to apply the A BCD law. We have 

[sin nd — sin (n — l)0]oi + d sin nd . . 

qi = ~ 7 — a : • (90) 

— t- sin nd-qi + I 1 — j J sin nd — sin (n — 1)6 

From (34) it follows that 

sin nd — sin (n — 1)6 = (d/2f) sin nd + sin 6 cos nd (91) 



which can be used together with (89) to rewrite (90) as 

f_l_ B sinflV" 9 fi-i^ 

92 *• sin0.e-*> + ^-^sinn<A 91 ^ ' 

After some further rearranging this can be written in the form 
1 



(92) 



1 I + _2 1 



qi q m q m f _qi q m q m f _ 



1 _L 

1 1 + 2_ 1 



exp (- 2jn6). (93) 



For the case where the q parameter of the injected beam does not 
differ too much from the parameter q m of a mode we can put 

A = (l/ ft ) - (l/q m ) (94) 
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and assume that A is small. Developing (92) in powers of A we obtain 

L _ I = a-^"" - jA 2 2^ (e 2i " B - e iinB ) 4- 0(A 3 ). (95) 
<?2 q m 2X 

If we neglect all but the first-order term in (95) and compare the real 
and imaginary parts, we arrive at approximate formulae* for the output 
parameters R 2 and w 2 : 



' r- i =(w + £) cos 2 nfl + * (-L - i) sin 2 «*, 
2/ \i?i 2// tt \wr w m 2 / 



[96) 



Bs ' 2/ 
— - — = - -(— + i ] sin 2nd + ( — , - —J cos 2 n0. 

Comparing these expressions with (33) we see that the beam radius 
w 2 varies in z direction with a period that is half the period with which 
a ray displacement varies. This fact has already been seen experimen- 
tally, 7 and has also been noted for other optical structures. 

The formulae (96) are valid for cases where the parameters Wi and 
Ri of the input beam do not differ much from the parameters of the mode 
of the lens sequence (i.e. for small A). For cases where this condition is 
not fulfilled we have to go back to (90). Using (72) we re-express the 
q parameters in terms of Wi , Ri , w 2 , and R 2 and compare the imaginary 
parts of l/q 2 as given by (90). After some algebra, where (91) is used 
to make simplifications, we obtain 



W 2 

w{ 2 



lb-^-(^(k+$h 2 «° (97) 



+ 2 



+ 

In this exact expression for w 2 we find the same periodicity in z direction 
as in (96). As n is varied w 2 goes through maximum values w max and 
minimum values w min . It is easy to show from (97) that 

u>mas w min = wj. 

Note that w mas and w min are the extrema of the envelope curve obtained 
for continuously variable n. The extremal values of w 2 actually occur at 
a lens only if the corresponding n is an integer. 

* In a recent publication by J. Hirano and Y. Fukatsu in Proc. IEEE, 52, Nov., 
1964, p. 1284, similar expressions were derived by means of a perturbation tech- 
nique in which the real beam parameters were used directly. 
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An exact expression for R% is obtained by comparing the real parts of 
l/g 2 in (90) in a similar way. 

4.7 Beam Transformation by a Lenslike Medium 

The passage of Gaussian beams through a lenslike medium as de- 
scribed by (42) has been discussed by several investigators. 19,28,29 ' 30,31 
We assume here for simplicity that n Q = 1, or a refractive index given by 

n = 1 - 2(r 2 /b 2 ). (98) 

It is easy to show " that for a Gaussian beam that is injected 

with a plane wave front and a beam radius w given by 

wo 2 = \b/2ir (99) 

the wave front remains plane, and the beam radius remains constant 
as the wave propagates. These light beams are called the modes of the 
lenslike medium. If the beam is injected with a beam radius W\ 9^ w , the 
wave front and the beam radius will change as a function of z. This 
problem has been treated by Tien et al. 28 with the help of a differential 
equation, and by Marcatili 29 who expanded the field distribution of the 
input beam in terms of the modes of the lenslike medium. We will show 
here that one can get the desired results in a rather simple fashion by 
employing the ABCD law (80). 

The elements of the ray matrix of a medium section of length z are 
given in (44). Using these together with (80) one computes for a beam 
with the complex parameter qi in the input plane a beam parameter 

z . b . z 
qi cos 2=- + -sin 2- 

q 2 = b 2 L (100) 

2.2 z 

— ffi r sin 2r + cos 2t 
b b 

in the output plane a distance z away from the input. Assuming an 
input beam with a plane wave front and a beam radius W\ we have 

-= -j— 9 - (101) 

qi irWi 2 

Inserting this and (99) in (100) we obtain 



. n z .Wo 2 n z 

1 ^ sm2 i; + 3—1 \ cos2 u 

1 _ __ _A W1 2 

§2 TTW0 2 n Z -Wo 2 • n Z 

COS 2r - J — sm 2r 



(102) 



OPTICAL MODES 470 

If we compare the imaginary parts in this expression we get 

wi = w? (cos 2 2- u + 2L sin 2 2?) (103) 

\ Wi* 0/ 

which agrees with the results of Refs. 28 and 29. A comparison of the 
real parts yields an expression for the curvature of the wave front. 

V. RESONATORS WITH INTERNAL OPTICAL ELEMENTS 

5.1 The Basic Resonator Parameters 

A resonator consisting of two spherical mirrors spaced at a distance 
d is shown in Fig. 12. Si and Z? 2 are the radii of curvature of the two 
mirrors, measured positive as shown in the figure. The mirror diameters 
or widths are 2a x and 2a 2 , respectively. The three basic parameters of 
such a resonator are 10 - 32,33 



am 

\d ' 



*-=0-«)' 



ft-£(l 



d 



(104) 
(105) 
(106) 



Within the Fresnol diffraction theory of optical resonators these three 




Fig. 12 — Empty spherical mirror resonator. 
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parameters determine completely the diffraction losses, the resonant 
frequencies, and the mode patterns of the resonator.' 

In the following we will show that resonators in which lenses or similar 
optical structures are inserted between the resonator mirrors are equiva- 
lent to an empty resonator of the type shown in Fig. 12. By equivalent 
resonators we mean here resonators with the same diffraction losses, 
the same mode patterns except for a scaling factor, and the same resonant 
conditions. To specify an empty resonator equivalent to a resonator 
with interna] optical elements we will compute its parameters N, Gi , 
and G 2 . 

5.2 Resonators with an Internal Lens 

A resonator with an internal lens is shown in Fig. 13. A lens of focal 
length / is spaced a distance di away from the left mirror and d% away 
from the right mirror. As before we call the radii of curvature of the 
two mirrors Ri and Rz , and their diameters 2a a and 2a 2 as shown. The 
internal lens is assumed to be so large that no additional aperture 
diffraction effects are introduced. 

Suppose now that we know the modes of the resonator. We can apply 
the imaging rules of Section II and choose the mirror surface of the right 
mirror, say, as reference surface. The image of the mode pattern on this 
mirror appears a distance 




Fig. 13 — Resonator with internal lens and equivalent empty resonator. 



OPTICAL MODES 481 

di = J^-. (107) 

d 2 - f 

away from the lens as shown in the figure. The field of the wave reflected 
from the mirror is zero outside the mirror aperture a 2 . The field of the 
corresponding image is therefore zero outside an aperture a 2 ' given by 
the magnification 

%= -%=\- d -}. (108) 

The image is a scaled reproduction of the mode pattern on the mirror 
which is exact in amplitude and phase if a spherical reference surface 
is chosen. The correct curvature of this surface is found in accordance 
with (6) and (8) by imaging the center of curvature of the mirror on the 
right. 

Consider now a mirror of diameter 2a 2 ' placed at the location of the 
image a distance 

d = ch - d 2 ' (109) 

away from the original left mirror as shown in the lower part of Fig. 13. 
The mirror curvature is chosen to be the same as the curvature of the 
reference surface for the image. This mirror may be called the- image 
mirror of the original mirror on the right. Apart from a phase difference 
of 2k (d 2 + d 2 ' ) it reflects a wave coming in from the left in exactly the 
same way as the original mirror combined with the Tens. The incoming 
wave produces the same (magnified) complex amplitude distribution 
or field pattern on the image mirror as on the original mirror on the 
right. The outgoing wave reflected by the image mirror has a field pattern 
at d 2 ' that is identical to the field pattern at d 2 of the outgoing wave 
reflected by the combination of the original mirror and the lens. The 
field patterns of the two outgoing waves are thus also identical in any 
other beam cross section, and in particular across the left mirror. There- 
fore the modes of the empty resonator formed by the image of the right 
mirror and the original left mirror as shown in the figure are equivalent 
to the original resonator with the internal lens. The mode patterns on 
the left mirror are identical in both cases, and the mode patterns of the 
corresponding mirrors on the right are scales of each other. The diffrac- 
tion losses of the two systems are also the same, and there is only a small 
difference in the corresponding resonant conditions due to the difference 
in phase shift of k (tf 2 + d 2 ) per transit. 

The basic parameters of the equivalent empty resonator are easily 
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obtained. According to (104) we have a Fresnel number of 

N = aia 2 '/\d (110) 

and with (105) 

Gl =%(l-±). (Ill) 

With (107), (108), and (109) these expressions can be written in terms 
of the dimensions of the resonator with the internal lens. One obtains 

N = S ^ ^x , (H2) 






X( di + d 2 - -V 

and 

d\ | , d 2 1/7,7 d\d 2 



a -5 1 -7-5V* + A -T ) ^ <U3) 



By interchanging subscripts one gets 



a 2 J , di 1/7,7 ^1^2 



< *-s( i -7-M* + *-7 ! ;t- (114) 

These three parameters determine the properties of the modes of the 
internal lens resonator. In the above expression one notes the appearance 
of the term 

d = di + di - (dA/f) (115) 

which one might call the effective distance between the mirrors. It is 
modified by the presence of the lens. 

In Refs. 4 and 32 approximate expressions are given for the resonant 
condition and the beam radii (spot size) of the fundamental mode at 
the mirrors of an empty resonator that is stable. Recall that for a stable 
resonator there holds 

g GiGt ^ 1. (116) 

We can apply these formulae to our equivalent resonator and obtain 
by imaging the corresponding expressions for the resonant wavelength 
X and the beam radii w t and w 2 on the mirrors of our resonator with an 
internal lens. Using the parameters discussed above we get 
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2(dl + d2) = ff +-(m + »+l) cos -1 VM2 (117) 

X IT 

where 9 is the longitudinal mode number, and m and n are the transverse 
mode numbers. The sign of the square root should be chosen equal to 
the sign of G { (or Q 2 ). For the beam radii we get 

wm = M?(i -G l G 2 )~\ (118) 



and 



©'• 



!* = ^(^?] . (119) 

The image mirror discussed above can also be obtained from the 
concept of a "thick mirror." A thick mirror 13 is a combination of a 
spherical mirror and a lens. The optical characteristics of this combina- 
tion are represented by a combined focal length and a principal plane. 
A mirror of this focal length located at the principal plane is equivalent 
to the thick mirror combination. This equivalent mirror is the same as 
our image mirror. 

The equivalence of the empty resonator and the internal lens resonator 
can also be shown by using the Fresnel diffraction formula in the manner 
of Appendix A. One obtains integral equations for the modes of an 
internal lens resonator. After performing the integration over the lens 
plane which involves infinite Fresnel integrals, the equivalence to the 
empty resonator is easily seen. 

For cases where the effective distance do as given by (115) is very 
small, ray angles of interest become rather large and the theory of 
Fresnel diffraction is no longer expected to apply. We have to exclude 
these cases from our considerations. 

Our discussion includes internal lens resonators with flat mirrors as 
shown in Fig. 14. The basic parameters of this resonator type can be 
obtained from (112), (113), and (114) by putting ffi = R* = co.Burch 
and Toraldo di Francia 8 have discussed the confocal system of this 
resonator family where G t = G 2 = 0. The transmission line dual of an 
internal lens resonator with flat mirrors is also shown in Fig. 14. It is a 
sequence of lenses and irises spaced as shown. In this sequence the lenses 
are large and the irises inserted between them control the modes of the 
system. For a symmetric system of this kind where di = d 2 = d and 
ai = a,2 = a the above expressions simplify, and we have 

ft = G, = 1 - (d/f), (120) 
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Fig. 14 — Internal lens resonator with flat mirrors and its transmission line 
dual, a sequence of lenses with irises placed between the lenses. 

and a Fresnel number of 



iV = 



(-»' 



\d 



(121) 



5.3 Resonators with an Internal Optical System 

As discussed before, the imaging rules of Section II apply not only 
to thin lenses but to any optical system that can be characterized by a 
focal length and by principal planes. The expressions derived above for 
internal lens resonators can therefore be applied also to spherical mirror 
resonators with an internal optical system. All one has to do is to inter- 
pret / as the focal length of the system and put 



ch = hi 



a\ = ho 



(122) 



where hi and hi measure the distances between the two principal planes 
and the corresponding mirrors. 

We can also characterize the internal optical system by its A BCD 
or ray matrix as in (9). Inserting (122) in (112), (113), and (114) we 
compare the resulting expressions with (19) through (22). We note 
immediately that the three basic resonator parameters can be written 
in terms of the elements of the ray matrix in the form 



M = qiCt2 
"" \B 



(123) 
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*-2("-$)- 



(124) 
(125) 



5.4 Internal Lenslike Medium — Guiding Medium with Apertures 

In this section we consider a spherical mirror resonator with a lenslike 
medium inserted between the resonator mirrors. The optical properties 
of a lenslike medium have been discussed in Sections 3.4 and 4.7. The 
refractive index of this medium changes with the square of the distance 
from the optic axis and is described by (98) if we assume n = 1. The 
degree of this index variation is measured by the parameter b. As shown 
in Fig. 15, we assume that the medium fills the space between the 
resonator mirrors which arc spaced at a distance I. The mirror diameters 
are 2a x and 2a 2 , respectively, and the corresponding radii of curvature 
are R\ and 7? 2 . 

The three basic resonator parameters which describe the modal prop- 
erties of this resonator with an internal lenslike medium are easily 
computed by using the results obtained before. The elements of the ray 
matrix for a medium section of length I are given in (44). Inserting 
these in (123), (124), and (125) we obtain for the Fresnel number of 
the system 



N = 



2aia 2 



\b sin 2 [ ' 

o 



(126) 



and for the G parameters 




r-K— - 



Fig. 15 — Resonator with an internal lenslike medium. 
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' <- 4-i sin2 l 



a 2 \ 

h = * (< 



■■ ' "4-^™n, 



i). 
0- 



(127) 
(128) 



A special case of the above system is shown in Fig. 16, where the 
mirrors are flat, i.e., R x = R 2 = <». The transmission line dual of this 
resonator is also shown in the figure. It is the interesting case of a lenslike 
medium or a gas lens with periodically spaced irises as shown. For irises 
of equal diameter (a\ = a> = a) the above expressions simplify, and 
we obtain for the Fresnel number of the system 

2a 2 

(129) 



N = 



and 



\b sin 2 i 
o 



Gi = Gi = cos 2 T . 
b 



(130) 



This system is confocal for I = (ir/4:)b. When the value of 2l/b ap- 
proaches a multiple of t, N gets very large and we have a case where the 
effective distance between the mirrors is very small [compare (115)]. 
As discussed before, the theory of Fresnel diffraction is no longer expected 
to apply under these circumstances. 
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Fig. 16 — Resonator with flat mirrors and an internal lenslike medium, and 
its transmission line dual, a gas lens with periodically spaced irises. 
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In high-gain lasers the parameter l/b can become rather large for 
certain frequencies. 21 For frequencies where the laser medium is focusing 
we have b 2 > 0, while for frequencies where the medium is defocusing 
b 2 < and b is imaginary. It is interesting to study the stability of a 
resonator with an internal lenslike medium allowing for positive and 
negative values of b 2 . For simplicity we assume that the radii of curvature 
of the two resonator mirrors are equal and put Ri = R-> = 2/. With 
(127) and (128) we obtain 



G,(? 2 = G 2 = (cos 2 l - - A sin 2 l - 

and write the stability condition (116) in the form 

-1 < G g 1. 



(131) 



(132) 



One can plot a stability diagram in which each resonator with given 
parameters l,f, and b is represented by a point. Such a diagram is shown 
in Fig. 17, where l/f is plotted as ordinate and l/b and jl/b are plotted 
as abscissae. Resonators with b 2 > are represented by points to the 




Fig. 17 — Stability diagram for a resonator with an internal lenslike medium. 
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right of the l/f axis, and resonators with b~ < are represented by points 
to the left. Points in the shaded regions correspond to unstable reso- 
nators, and resonators represented by points in the unshaded regions 
are stable. The boundaries between the stable and unstable regions 
follow from (131) and (132). They are described by the equations 

j- = v ^, v integer, (133) 

j=4- & cot^, (134) 

l -= -4J-tenL (135) 

* bo 

For b 2 < 0, where b is imaginary, the trigonometric functions of (134) 
and (135) become hyperbolic functions as in (48) and (49). For b > 
one gets periodically stable and unstable regions as l/b is increased. 

We have not discussed in detail cases where the lenslike medium 
occupies only a part of the space between the resonator mirrors. How- 
ever, one can compute easily the basic parameters for resonators of this 
kind with the help of the matrix elements of (44), and the formulae 
(123), (124), and (125). 

5.5 Resonators with One Very Large Mirror 

Let us return to the case of an empty resonator. In some practical 
arrangements the diameter of one of the two mirrors, say 2a.2 , is so 
large that diffraction by its aperture can be neglected. The resonator 
modes are then more or less controlled by the aperture a x of the other 
mirror. This statement is not true for resonators of the degenerate 
confocal type where the diffraction losses at each mirror are equal for 
any aperture ratio a 2 /ai . We exclude resonators of this type from our 
present discussion. 

The properties of the resonator modes are generally determined by 
the three basic parameters given in (104), (105), and (106). But for 
an infinitely large a<t the Fresnel number N and the parameter Gi be- 
come infinitely large, and G\ = 0. The resonator parameters are now 
quite meaningless. It is, however, possible to construct an equivalent 
resonator with parameters of finite value, as we will show below. 

Consider Fig. 18. An empty resonator with one mirror of large diameter 
is shown schematically at the top. Below it we have drawn its transmis- 
sion line dual. It consists of a sequence of lenses where an apertured lens 
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Fig. 18 — Empty resonator with one very large mirror, its transmission line 
dual, and its equivalent internal lens resonator. 

follows an unapertured lens of large diameter. But this transmission 
line is also the dual of the resonator shown at the bottom of tiie figure. 
This is a resonator formed by apertured mirrors of finite diameter 2a x 
with an internal lens of focal length / = R 2 /2. The lens is unapertured. 
Internal lens resonators of this type have been considered before. We 
can compute the Fresnel number of this system from (112) and obtain 



N = 



ar 



2Xd[l - ■£■ 



(-4)' 



(136) 



Equations (113) and (114) are used to calculate the G parameters with 
the result 



G x = G 2 = 1 - 



Ki + k-ik)- (137) 
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These parameters determine the properties of the modes of the internal 
lens resonator shown in Fig. 18. The mode patterns at the apertured 
mirrors of this resonator are, of course, equal to the mode pattern at the 
apertured mirror of the empty resonator. The one-trip diffraction loss 
of a mode of the internal lens resonator is equal to the return-trip diffrac- 
tion loss of an empty resonator mode, as there are no diffraction losses 
at the infinitely large mirror. 

For the special case where the large mirror is flat CR 2 = °° ) the above 
discussed equivalences are well known. They follow from symmetry 
considerations. 
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APPENDIX A 

Imaging for Fresnel Diffraction 

The purpose of this appendix is to show how the imaging relation (6 ) 
of the main text is derived within the formalism of scalar Fresnel diffrac- 
tion theory. Assume a light wave traveling in z direction and refer to 
Fig. 19. Call the object field Ei(xi , yi) and the image field Es(x} ,yz). 
The distances d\ and di between the lens and the object and image 
planes, respectively, are related by 

(1/di) + (l/dt) = 1/f. (138) 



OBJECT PLANE 
'I. 






X 

+■ 

L 



IMAGE PLANE 
X 2 

* 



f_ 

Ei tamo Eodoiyo) Eo(ioiyo) E 2 (i 2 ,y 2 ) 

Fig. 19 — Dimensions of interest for Fresnel diffraction theory of imaging. 
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The field immediately to the left of the lens is denoted E u (x { , , y ) and 
the field to the right of* the lens is E (x a ,Uo)- According to (1) of the 
main text we have for a large, ideal lens 

E = E»exp(-jk X ±±^ (139) 

where k = 2w/\ is the propagation constant in the medium. With the 
help of the Fresnel diffraction formula the fields En and E- 2 can be ex- 
pressed as 

E = JL f dxxdyxEi exp (-jfepi) (140) 

27T«i >>A 1 

E 2 = J^- [ dxod.yoEn exp (-jkp 2 ) (141) 

ZtTUi J- oo 



and 



where 



and 



Pl = dl + 2^ (X1 " Xo)t + 2^ iVl " Vo)2 iym 



p 2 = d, 4- ± fa - .to) 2 4- ±r (1J2 - 2/o) 2 . (143) 

The integration in (140) is performed over the aperture area A x of the 
object field, and the integration limits in (141) are extended to infinity 
with the assumption that the lens is so large that no additional aperture 
diffraction effects are introduced. 

Combining (139), (140), and (141) we obtain by interchanging the 
order of integration 

k 2 f f + °° 

Ei= - . ', , / dxidyxEi / dxtfhji) 

•exp 

where 

r 2 = .ro 2 + ?/0 2 . (145) 

Now the expressions (142) and (143) for pi and p 2 are inserted. One 
finds that in the exponential the terms proportional to r cancel because 
of (138). The integration with respect to x and ?/o can be performed by 
noting that 



-jk ( Pl + P2 4- ? ^jj 



(144) 
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£" <te exp [;,,,„ (| + g] = W (* [| + |]) ( 146) 
where 5 is the Dirac delta function. 34 With this (144) becomes 

• Jf dgrfftg, exp (-# jkj (147) 

■•('[i+^ft+SD- 

This simplifies immediately with the help of the formalism of the delta 
function to 

■-p[-*(* + * + a( 1 + S)]- 

Multiplying (138) by d\/a\ one finds that 

3-6+D-7J (149) 

di \ d 2 / f a* 
which is used to write (148) in the form of (6) of the main text. 

APPENDIX B 

Principle of Equal Optical Path Leading to Additional Phase Shift in 
Image Plane 

The process of imaging the field distribution in the object plane into 
the image plane can be understood in terms of the rays leaving each 
point (say Pi) in the object plane at various angles as shown in Fig. 20. 
All rays originating from Pi are collected at a corresponding point P 2 
in the image plane. A form of the principle of equal optical path 36 says 
that the optical path lengths from Pi to P 2 are the same for all rays 
regardless of initial slope. 

To obtain an image which is an exact reproduction of the original 
amplitude and phase distribution it would be necessary for the various 
optical paths which connect corresponding points, say Pi and P 2 or 
Qi and Q 2 , to be equally long for all points regardless of their distance 
from the optic axis. That these path lengths are not the same for all 
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OBJECT PLANE 



IMAGE PLANE 




Fig. 20 — Rays emerging from a point of the object collected at the image. 

points but increase with increasing distance between the imaged point 
and the optic axis can be seen from the simple example of an ideal plane 
wave coming in from the left. This case furnishes an expression for the- 
path length difference as a function of the distance between, the imaged 
point and the axis. As the path length is independent of the field distribu- 
tion imaged, this expression is valid for the general case. To derive it 
we recall that an ideal plane wave is transformed by an ideal lens into 
an ideal spherical wave with the focal point of the lens as its center. 
The rays connecting points which lie on corresponding wave fronts are 
equally long for all points on the wave front. 35 Therefore all path lengths 
measured from the object plane to the spherical wave front which touches 
the image plane are equal. Paraxial rays (which are practically parallel 
to the optic axis) need an additional length equal to r 2 2 /2f to reach a 
point (P 2 ) in the image plane which is a distance r 2 away from the axis. 
This additional ray length accounts for the additional phase shift given 
in (6) in the main text. 



REFERENCES 

1. Fox, A. G., and Li, Tingye, Resonant Modes in a Maser Interferometer, 

B.S.T.J., 40, March, 1961, p. 453; Proc. IRE, 48, I960, p. 1904. 

2. Goubau, G., and Schwering, F., On the Guided Propagation of Electromag- 

netic Wave Beams, Trans. IRE, AP-9, May, 1961, p. 248. 

3. Boyd, G. D., and Gordon, J. P., Confocal Multimode Resonator for Milli- 

meter Through Optical Wavelength Masers, B. S.T.J. , 40, March, 1961, p. 
489. 

4. Boyd, G. D., and Kogelnik, H., Generalized Confocal Resonator Theory, 

B.S.T.J., 41, July, 1962, p. 1347. 

5. Fox, A. G., and Li, Tingye, Modes in a Maser Interferometer with Curved 

and Tilted Mirrors, Proc. IEEE, 51, Jan., 1963, p. 80. 

6. Goubau, G., Optical Relations for Coherent Wave Beams, in Electromagnetic 

Theory and Antennas, ed. Jordan, E. C, Pergamon Press, 1963, p. 907. 

7. Kogelnik, H., Matching of Optical Modes, B. S.T.J. , 42, Jan., 1964, p. 334. 



494 THE BELL SYSTEM TECHNICAL JOURNAL, MARCH 1965 

8. Burch, J. M., Design of Resonators, p. 1187 in Quantum Electronics III, 

ed. Grivet, P., and Bloembergen, N., Columbia University Press, New 
York, 1964; Toraldo di Francia, G., On the Theory of Optical Resonators, 
Proc. Symp. Optical Masers, New York, 19(53, p. 157, Polytechnic Press, 
Brooklyn, New York. 

9. Kogelnik, H., and Yariv, A., Considerations of Noise and Schemes for Its 

Reduction in Laser Amplifiers, Proc. IEEE, 52, Feb., 1964, p. 165. 

10. Gloge, D., Analysis of Fabry -Perot Laser-Resonators by Means of Scatter- 

ing Matrices, Arch. El. U., 18, March, 1964, p. 197. 

11. Collins, S. A., Analysis of Optical Resonators Involving Focusing Elements, 

Appl. Opt., 3, Nov., 1964, p. 1263. 

12. Li, T., Dual Forms of the Gaussian Beam Chart, Appl. Opt., S, Nov., 1964, 

p. 1315. 

13. Jenkins, F. A., and White, H. E., Fundamentals of Optics, 3rd ed., McGraw- 

Hill Book Company, New York, 1957, p. 89. 

14. Born, M., and Wolf, E., Principles of Optics, Pergamon Press, New York, 

1959, p. 420. 

15. Pierce, J. R., Theory and Design of Electron Beams, D. Van Nostrand Com- 

pany, Inc., New York, 1954, p. 194. 

16. Herriott, D. R., Kogelnik, H., and Kompfner, R., Off-Axis Paths in Spherical 

Mirror Interferometers, Appl. Opt., 3, April, 1964, p. 523. 

17. Bertolotti, M., Matrix Representation of Geometrical Properties of Laser 

Cavities, Nuovo Cimento, June, 1964, p. 1242; O'Neill, E. L., Introduction 
to Statistical Optics, Addison-Wesley Publ. Co., Reading, Mass., 1963; 
Brouwer, W., Matrix Methods in Optical Instrument Design, Benjamin, 
New York, 1964. 

18. Richards, P. I., Manual of Mathematical Physics, Pergamon Press, London, 

1959, p. 312. . . , 

19. Tonks, L., Filamentary Standing-Wave Pattern in a Solid-State Maser, J. 

Appl. Phys., June, 1962, p. 1980. 

20. Gudzenko. L. I., Concentration of a Light Wave in a Weakly Inhomogeneous 

Dielectric, J. Exp. Theor. Phys. (USSR), 44, April, 1963, p. 1298; Sov. 
Phys. JETP, 17, Oct., 1963, p. 875. 

21. Kompfner, R., Talk at the Electron Device Research Conference, Salt Lake 

City, Utah, 1963. 

22. Berreman, D. W., A Lens or Light Guide Using Convcctively Distorted 

Thermal Gradients in Gases, B.S.T.J., 43, July, 1964, p. 1469; A Gas Lens 
Using Unlike, Counter-Flowing Gases, B. S.T.J. , 43, July, 1964, p. 1476. 

23. Marcuse, D., and Miller, S. E., Analysis of a Tubular Gas Lens, B. S.T.J. , 

43, July, 1964, p. 1759. 

24. Beck, A. C, Thermal Gas Lens Measurements, B.S.T.J., 43, July, 1964, p. 

1818; Gas Mixture Lens Measurements, B.S.T.J., 43, July, 1964, p. 1821. 

25. Born and Wolf, op. cit., p. 121. _ 

26. Kogelnik, H., Modes in Optical Resonators; in Advances in Lasers, ed. Levine, 

A. K., Dekker Publishers, New York. 

27. Yariv, A., and Gordon, J. P., The Laser, Proc. IEEE, SI, Jan., 1963, p. 4. 

28. Tien, P. K., Gordon, J. P., and Whinnery, J. R., Focusing of a Light Beam 

of Hermite-Gaussian Distribution in Continuous and Periodic Lenslike 
Media, to be published in Proc. IEEE. 

29. Marcatili, E. A. J., Modes in a Sequence of Thick Astigmatic Lens-Like 

Focusers, B.S.T.J., 43, November, 1964, p. 2887. 

30. Gloge, D., Focusing of Coherent Light Havs in a Space Dependent Dielectric; 

Arch. El. V., 18, July, 1964, p. 451. 

31. Pierce, J. R., unpublished work. 

32. Gordon, J. P., and Kogelnik, H., Equivalence Relations among Spherical 

Mirror Optical Resonators; B.S.T.J., 43, November, 1964, p. 2873. 

33. Streifer, W., and Gamo, H., On the Schmidt Expansion for Optical Resonator 

Modes, Proc. Symp. on Quasi-Optics, Polytechnic Inst, of Brooklyn, June, 
1964. 

34. Born and Wolf, op. cit., p. 752. 

35. Born and Wolf, op. cit., p. 126. 



